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Abstract 

We examine various direct and indirect experimental constraints on the nearly minimal super- 

Oh; 

symmetric standard model (nMSSM) and obtain the following observations: (i) Current experi- 
ments stringently constrain the parameter space, setting a range of 1 ~ 37 GeV for the lightest 
neutralino Xi> 30 ~ 140 GeV (1 ~ 250 GeV) for the lightest CP-even (CP-odd) Higgs boson, and 
1.5 ~ 10 for tan /3; (ii) To account for the dark matter relic density, besides the s-channel exchange 



> 

^ \ of a Z-boson, the s-channel exchange of a light A\ (the lightest CP-odd Higgs boson) can also play 



o 

O 



an important role in LSP annihilation. Compared with the Z-exchange annihilation channel, the 
A\ exchange channel is more favored by muon g — 2 data and allows much broader regions for 
the parameters; (iii) In a large part of the allowed parameter space the SM-like Higgs boson may 
dominantly decay to XiXi or ^1^1 an d the conventional visible decays (e.g. into bottom quarks) 
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are severely suppressed. 



PACS numbers: 14.80. Cp,12.60.Pr,ll. 30. Qc 



1 



Because the minimal super symmetric standard model (MSSM) suffers from the /i- 
problem, some non-minimal supersymmetric models have recently been intensively stud- 
ied, amongwhich an attractive one is the nearly minimal supersymmetric standard model 
(nMSSM) Jl[ Q]. This model extends the MSSM by one singlet superfield S with the super- 
potential [3[] 

W = W M ssm + teiiKHiS + &M*S, (1) 

where Wmssm is the superpotential of the MSSM without the /x-term, the second term 
on the right side is the interaction of the singlet S with the Higgs doublets H u and Hd, 
and the last term is the tadpole term. This superpotential differs from that of the next- 
to-minimal supersymmetric standard model (NMSSM) 4| in that the tadpole term of the 
nMSSM replaces the trilinear sing let term kS 3 of the NMSSM. Due to this tadpole term, W 
has no discrete symmetry and so the nMSSM is free of the domain wall problem suffered by 
the NMSSM. The tadpole term also gives rise to a vacuum expectation value (vev) for the 
singlet, controlled by £,fM%. Though it is SUSY-preserving, this tadpole term can naturally 
be of the SUSY breaking scale - e.g., in the N — 1 supergravity model with a discrete 



R-symmetry [20(] such a tadpole term is generated at a high loop level and thus is naturally 
small [3]. A nonzero singlet vev at the SUSY breaking scale generates an effective /x term 
from the \EijH l u H d S term with the desired order of magnitude, solving the /x problem of the 
MSSM. These theoretical virtues motivate further phenomenological study of the nMSSM. 
Because of the absence of the trilinear singlet term, the spectrum and phenomenology of the 
nMSSM can be quite different from those of the NMSSM. 

With the running of the LHC, all low energy supersymmetric models will soon be put 
to the test. To explore these models at the LHC, it is very important to determine the 
parameter space allowed by current experiments. In this work we comprehensively examine 
experimental constraints on this model from the direct experimental searches for Higgs 
bosons and sparticles, the precision electroweak measurements at LEP/SLD, the cosmic 
dark matter relic density from WMAP, and the muon anomalous magnetic moment. We 
also consider the theoretical constraints from the stability of the Higgs potential and the 
perturbativity of the theory up to the grand unification scale. After analyzing the allowed 
parameter space, we discuss some phenomenology of this model. 

We start our analysis by recapitulating the basics of the nMSSM. With the superpotential 



in Eq. ([I]), the corresponding soft-breaking terms are given by 

V soit = Vmssm + ml\H d \ 2 + ml\H u \ 2 + m 2 s \S\ 2 

+(XA X £ ij HiH 3 d S + £ S M*S + h.c.) (2) 

where Vmssm contains the soft breaking terms for gauginos and sfermions in the MSSM, 
and m U: d,Si A\ and £sM^ are soft breaking parameters. Noting that the tadpole terms do 
not induce any interactions, one can conclude that, except for the tree- level Higgs boson 
masses and the minimization conditions, the theory is the same as the well-known NMSSM 
with the trilinear singlet term set to zero lal. The nMSSM predicts three CP-even and two 

n n 

CP-odd neutral Higgs bosons as well as five neutralinos [1|, |5j. The mass of the lightest 
neutralino (assumed to be the LSP) arises from the mixing of the singlino with higgsinos 
and is given by [6j] 

2fi\ 2 (vl + v 2 d ) tan/3 



y > 2/i 2 + \ 2 {v 2 u + v 2 )tcm 2 (3 + 1 



(3) 



where tan/3 = v u /v d , fi = A(s), and v u , v d and (s) are the vevs of the Higgs fields H u , H d 
and S, respectively. 

In our calculations we extend the packages NMSSMTools [7] and micrOMEGAs js] to 
the nMSSM. We use the modified NMSSMTools to calculate the Higgs boson masses and 
their decays including all known radiative corrections. We use the modified micrOMEGAs 
to calculate the dark matter relic density. The parameters relevant to our analysis are A, 
A\, tan/3, uia = 2{^iA\ + \£pM 2 )/ sin2/3, /1 = X(s), rhs, the gaugino masses Mi and M 2 , 
and the soft SUSY breaking parameters in the squark/slepton sectors. We assume all these 
parameters to be real and specify their values at the weak scale. 

Since we are only interested in the properties of the Higgs bosons and neutralinos and 
these properties are affected little by the soft parameters in the squark/slepton sectors, we 
specify these parameters before our scan. Noting that a heavy stop (t) is helpful for the Higgs 
sector to evade the LEP constraints and a light smuon (jl) is needed in the nMSSM to explain 
the muon anomalous magnetic moment a M , we assume all the soft breaking parameters (soft 
masses and trilinear A-parameters) to be 1 TeV for the (t, b) sector and 100 GeV for the 
(z/^,/2) sector. We will briefly discuss the results with a lower mi and those with different 
rrifi when the a M constraint is switched on. For the other soft breaking parameters in the 
squark/slepton sector, we uniformly set them to be 1 TeV since the considered constraints 



are not sensitive to them. Moreover, we assume the grand unification relation for the gaugino 
masses, M\ = (5g' 2 /3g 2 )M2. With the above assumptions, we scan over the remaining seven 
parameters in the following ranges: 0.1 < A < 0.7, 1 < tan/5 < 10, —1 TeV < A\ < 1 TeV, 
50 GeV <m A ,fi,M 2 <l TeV and < m s < 200 GeV. Note that in the nMSSM tan f3 > 10 
is not allowed by the dark matter relic density because m^o is suppressed by large tan (3 (see 
Eq. ([3])) and a light Xi is difficult to annihilate sufficiently. 

In our scan we consider the following constraints: (1) The dark matter relic density, 
0.0945 < Vth 2 < 0.1287 [9J. We require Xi to account for this density. (2) The constraint, 
a e * p - a s ^ M = (29.5 ± 8.8) x 10~ 10 [10]. We require the nMSSM contribution to explain the 
deviation at the 2a level. (3) The LEP bound on invisible Z decay, T(Z — > XiXi) < 1-76 



MeV 



llj (we also apply this bound to the decay Z — > h\A\ with h\ being the lightest CP- 



even Higgs boson); the LEP-II upper bound on a(e + e~ — > XjXj)i which is 10~ 2 pb for i — 1, 
j > 1 and 10 _1 pb for i, j > 1 (summed over i and j) [12 ; and the lower mass bounds on 
sparticles from direct searches at LEP and the Tevatron [ll| . (4) Constraints from the direct 
search for Higgs bosons at LEP-II, which limit all possible channels for the production of 
Higgs bosons. (5) Constraints from precision electroweak observables such as pi ep t, sm2 ^ e //> 
my/ and R^. (6) The perturbativity of the nMSSM up to the grand unification scale and 
the stability of the Higgs potential which requires that the physical vacuum of the Higgs 
potential is the true minimum of the potential. 

The above constraints have been encoded in NMSSMTools Q], except for (5). In Ref. [3] 
we extended the code by adding (5); here we extend all these constraints to the nMSSM 
scenario. Since the hadronic contribution to a^ M remains under discussion [_lo]|, in the 
following we will present results both with and without considering the constraint. 



Our scan sample is 2.5 billion random points in the parameter space given above. With 
all the constraints except a M , only about 6 thousand points survive. This is mainly because 
the dark matter relic density stringently constrains the mass and couplings of the LSP, and 
consequently, only a small portion of the parameter space is allowed (as shown explicitly in 
Figs. 1-4 below). Among the surviving points, about 60% are characterized by > irtz 
and 30 GeV < m^o < 37 GeV, in which Xi mainly annihilates through Z-boson exchange to 
give the required dark matter relic density jlj . For most of the other points, both x\ an d A\ 
are predominantly singlet-like with roughly 2m^o ~ m^, in which the exchange of a light 
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Ai is the main annihilation channel of Xi (this annihilation channel, similar to the case of 
Table III in Q for the NMSSM, has not been studied for the nMSSM in the literature). 

If we switch on the a M constraint, the surviving points are reduced to 1.4 thousand, among 
which about 60% (70%) satisfy < 60 GeV (m Al < m z ). This indicates that the current 
a M measurement can constrain the model stringently, and that it favors LSP annihilation 
through exchange of a light A\ rather than a Z-boson. Note that in the above we fixed 



VTL; 



100 GeV. If we increase run, even fewer scan points will survive. For example, raising 



to 150 GeV eliminates about half the remaining points with < mz and nearly all the 
points with rriA 1 > %■ For m^ = 200 GeV, none of our scan points survive the a M constraint, 
which implies that smuons must be lighter than about 200 GeV. This conclusion is unique to 
the nMSSM. The underlying reason is that in SUSY models, the leading charging /neutralino 
contribution to 5a M is proportional to tan(3/m^ with N > 2 jl^ . [la ]. In the MSSM or 
NMSSM, tan (3 can be quite large [l3| and thus m^ is not stringently constrained by a M ; but 
in the nMSSM, tan/3 is bounded from above (^ 10) by the dark matter relic density and 
hence must be light. 
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FIG. 1: Parameter scan points allowed by current experiments. The dark bullets (light triangles) 
correspond to m^ a > mz (jriAi < m z), with Xi mainly annihilating through exchanging a Z-boson 
(a light A±) to give the required dark matter relic density. 
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FIG. 2: Same as Fig. 1, but showing dependence on fi, m^o, and m~+ 



In Figs. 1 and 2 we display the surviving points as functions of various parameters. As we 
pointed out before, the strongest constraint come from the dark matter relic density. These 
figures show that if only the Z exchange is responsible for the density (dark bullets), the 
parameters A, tan /3, fi, m^o 2 and m-+ are all constrained in quite narrow ranges; however, 
if the light A\ exchange is considered (light triangles), the allowed parameter ranges are 
significantly more spread out. 

Note that we checked by using the modified NMSSMTools that b-s transitions such as 
b — > sj and B®-B® mixing do not impose any meaningful constraints on the surviving 
samples in the case of no squark flavor mixings. We also checked that for samples with 
rriA 1 ^ 10 GeV, the branching ratio for T — ► ^Ai is less than its experimental upper bound 
of 1 x 10" 4 |ig|. 

The following additional comments are in order, (i) We call Xi singlino-like when the 
coefficient of the singlino component in x\ is larger than 1 / v2 (so its square is larger than 
1/2). In general the higgsino or gaugino components in Xi are n °t negligible even when 
Xi is singlino-like. In fact, it is the higgsino components in Xi that are mainly responsible 
for the LSP annihilation coupling discussed above. For annihilation via Z exchange, the 
typical coefficient of the if u -type higgsino component in Xi is 0.4, and we checked that as 
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the higgsino components increase, the coupling of XiXi^ increases and consequently the 
relic density drops [1]. For annihilation via A\ exchange, the typical coefficient for the 
higgsino component in Xi is 0.2 and for the doublet-Higgs component in A\ is 0.15. Since 
the couplings of an d Aiff (where / is a light fermion) are suppressed, annihilation 

via A\ exchange is too weak to produce the required relic density except in the funnel 
region 2m^o ~ jl^]. (ii) The requirement % < 200 GeV in our scan is taken from 
Ref. [l| which studied electroweak baryogenesis. We checked that a larger range of rhs does 
not change our conclusions qualitatively; it only increases the number of surviving samples 
with m J 4 1 > mz and raises the upper bound of m^ 1 . We did not impose the requirement 
of successful electroweak baryogenesis in our scan since there exist other ways to explain 
the origin of the matter-antimatter asymmetry in the Universe (iii) A heavy t is not 
necessary in the nMSSM for the SM-like Higgs boson to evade the LEP bound since the 
Higgs boson mass can be enhanced at tree level by a contribution from A [4]. We checked 
that a lighter t in our scan does not change our conclusions; it only decreases the number 
of the surviving points. 



As shown above, if we consider all the constraints including a M with = 100 GeV, the 
mass spectrum is limited to the following ranges: lGeV ^ ^ 250GeV, 30GeV ^ m hl ^ 
140GeV, 70GeV ^ rrih ^ 145GeV {h is the lightest doublet-dominant CP-even Higgs boson, 
usually called the SM-like Higgs boson; in some cases h can be hi), lGeV ^ m^o ^ 37GeV, 
50GeV ^ m^o ^ 300GeV and 105GeV ^ m~+ ^ 400GeV. In such allowed mass ranges, the 
phenomenology of the Higgs bosons and sparticles may be quite peculiar and different from 
the MSSM. A comprehensive study of the phenomenology of this model at the Tevatron and 
the LHC is beyond the scope of this paper; instead we present the following brief discussion. 

First, consider the Higgs bosons. The dominant decay of A\ may be either XiXi or bb for 
rriA 1 > 2mb, and we checked that the AiXiXi interaction is mainly induced by the higgsino 
components of Xi and/or the doublet components of Ai [7|. The dominant decay mode of 
the SM-like Higgs boson h can be any of the following: h — * XiXi^ X1X21 A1A1, hxh\. In Fig. 3 
we show the branching ratios of h and Ai decays into the LSP and bb. For ^ 60 GeV, 
both h and A% can decay predominantly into LSP pairs and the decay into bb is strongly 
suppressed. Since the branching ratio of h — > bb is suppressed below 10% in most of the 
allowed parameter space due to the presence of new competing decay modes, conventional 
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searches for a light SM-like CP-even Higgs boson at the LHC will be quite hopeless. For 
example, our results indicate that for about 60% of the surviving points, the final decay 



products of h are two or 
is a good search channel 



our LSPs. For these points, weak boson fusion with h — > invisible 
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Our results also indicate that for about 16% of the surviving 



points, h decays predominantly to 46. In this case, W/Zh production may be a good channel 



to detect h 



18|. 
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FIG. 3: Same as Fig. 1, but for the branching ratios of h and A\ versus m J 4 1 with all constraints 
imposed including cu. The effect of the constraint is to reduce the number of surviving points, 
as shown in Figs. 1 and 2. 

Second, consider smuon production at the LHC. If we take the constraint seriously, the 
smuon should be lighter than 200 GeV, which implies that it would be copiously produced 
either directly or from cascade decays of other sparticles at the LHC, and should be visible 
at the LHC. We checked that due to the non-negligible gaugino component of xl, the decay 
width of Jx — > fiXi is about several MeV, so the decay length of smuon is not macroscopic 
(for a heavy charged particle with macroscopic decay length, its signals at the LHC may be 
quite special [l^). 

Third, consider the next-to-lightest neutralino a ^ the LHC. As shown in Fig. 4, x\ i s 
bino-like for m^o ^ 150 GeV, and may be higgsino-like for m^o > 150 GeV. We checked that 
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FIG. 4: Same as Fig. 1, but showing dependence on Mi, m^o and I-/V21I with all constraints imposed 
including a M . N21 denotes the coefficient of the bino component in yQ. 

the main decay products of X2 can °e Xi^i, Xi^ij Xi^, or \\Z and there is a large portion 
of the surviving samples in which \% — * x\^i{hi-,h) — > XiXiXi- If X2 ^ s the NLSP and it 
mainly decays into 3x°, it will be copiously produced from the cascade decay of squarks and 
gluinos 2|, and can be easily mistaken as the LSP. From Fig.4 one can also learn that the 
gaugino mass Mi is bounded from 50 GeV to 500 GeV. This implies by the gaugino mass 
unification relation that the gluino mass varies from about 300 GeV to 3 TeV, which could 
be accessible at the LHC. 

In summary, we examined the current experimental constraints on the nMSSM. We found 
that the parameter space of this model is stringently constrained by current experiments, and 
in the allowed parameter space the phenomenology of this model may be quite peculiar. Such 
tightly constrained parameter space could make this model readily tested (either verified or 
excluded) at the LHC. In addition, since in this model the dark matter particle is constrained 
in a narrow mass range, the astrophysical dark matter experiments may also be able to cast 
some light on this model. 
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